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Abstract. Multiplicative cascades have been introduced in turbulence to generate random or deterministic 
fields having intermittent values and long-range power-law correlations. Generally this is done using discrete 
construction rules leading to discrete cascades. Here a causal log-normal stochastic process is introduced; 
its multifractal properties are demonstrated together with other properties such as the composition rule 
for scale dependence and stochastic differential equations for time and scale evolutions. This multifractal 
stochastic process is continuous in scale ratio and in time. It has a simple generating equation and can be 
used to generate sequentially time series of any length. 



PACS. 47.27. Jv High-Reynolds-number turbulence 
47.27.Eq Turbulence simulation and modeling 



1 Introduction 

Since their introduction in the fields of turbulence, geo- 
physics and chaos theory, multifractals have been widely 
used in many fields of applied sciences, including turbu- 
lence, geophysics, precipitations, oceanography, high en- 
ergy physics, astronomy, biology and finance. The main 
characteristics of multifractal models or data is to pos- 
sess high variability on a wide range of time or space 
scales, associated to intermittent fluctuations and long- 
range power-law correlations. In this framework, many 
studies have been devoted to data analysis and discrete 
mathematical model constructions. More precisely, there 
is a dichotomy in the multifractal literature between, on 
the one hand, experimental results showing usually no 
characteristic scales in a scaling range between an outer 
and an inner scale, and on the other hand, discrete cas- 
cade models that are not stationary and are scaling only 
for discrete scales ratios. This is why continuous multi- 
fractal models are needed, and among them, multifrac- 
tal stochastic processes: processes possessing multifractal 
statistics while depending on a continuous time parame- 
ter. Here we consider one such process, belonging to the 
generatic family of lognormal multifractals, while depend- 
ing on a continuous time parameter. This process has a 
simple analytical expression, and is causal and ID. We 
will discuss some of its properties, and confirm them us- 
ing data analysis of a numerical simulation. 

Below we first recall the basic framework of discrete 
cascades: their statistics and correlation properties. These 
properties are given as the definition of a multifractal cas- 
cade process, and are recovered for continuous causal pro- 
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cesses. In section 2, we discuss the continuous scale limit 
of discrete cascades, and provide explicitly a basic causal 
equation generating a multifractal stochastic process for 
log-normal multifractals. We explore their statistical prop- 
erties, including correlations for the process and its loga- 
rithm, and provide differential forms in time and in scale 
space. A link with the Langevin equation for Markov cas- 
cades is presented. Section 3 presents a simulation of a 
discrete cascade and a time simulation for a multifractal 
continuous process, together with a data analysis which 
confirms the theoretical results for both simulations. Sec- 
tion 4 is devoted to concluding remarks. Some useful re- 
sults in probability theory are recalled in the Appendices. 



1.1 The generic example of multiplicative cascades in 
turbulence 

We briefly recall here the classical framework of Richard- 
son turbulent cascades, which can be considered as the 
generic example of cascades. For high Reynolds number 
turbulence, dissipative terms are negligible compared to 
the inertial terms of the Navier-Stokes equations, so that 
for the inertial range, there is a local energy transfer from 
large to small scales. This is the picture of Richardson's 
1922 energy cascade ,1., which was later formalized by 
Kolmogorov in 1941 [2], corresponding to the Komogorov- 
Obukhov -5/3 Fourier power spectrum for velocity fluc- 
tuations in the inertial range. This description has been 
found experimentally valid for mean fluctuations, and this 
general cascade framework for the inertial range is now 
well accepted. Nevertheless it was not able to take into 
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Fig. 1. Illustration of the discrete cascade process. Each step 
is associated to a scale ratio of 2. After n steps, the total scale 
ratio is 2 n . 

account intermittent fluctuations of the small-scale dissi- 
pation field, as experimentally observed a few years later 
[3] . This intermittency lead Kolmogorov and Obukhov to 
propose new models in 1962 assuming log-normal 

statistics for the velocity increments and for the dissipa- 
tion field e, with furthermore the following asymptotic hy- 
pothesis for the variance of loge: 

^o ge ~^ + MlogA (1) 

where A and \i > are constants (fi is often called the 
intermittency exponent), and A = Ljl is the scale ratio 
of the outer scale L to the scale of interest I (in the same 
paper this property of the variance was also assumed for 
the averaged dissipation e r over a sphere of radius r) . Kol- 
mogorov's proposals in 1962 were given without any phys- 
ical or modelling justification; it was rather ad hoc. Soon 
after this, several experimental studies have shown that 
the dissipation field possesses long-range power-law corre- 
lations |rj][7j . This was a very specific property which was 
not included in the proposal of Kolmogorov. This lead 
Yaglom to propose in 1966 his explicit cascade model, 
which was able to reproduce experimental facts - small- 
scale intermittency and power-law correlations - and the 
log-normal hypothesis of Kolmogorov. This multiplicative 
cascade model had hence experimental grounds, and was 
proposing an unification of different results and proposals. 
It is discussed below. 



1.2 Discrete multiplicative cascades and their 
properties 

We present now discrete multiplicative cascades and their 
statistical properties. This is basically a discrete model 
in scale, but it can be densified. The terms "discrete in 
scale" refer to the fact that the scale ratio from mother 
to daughter structures is an integer number strictly larger 
than 1. This model is multiplicative, and embedded in 
a recursive manner. The multiplicative hypothesis gener- 
ates large fluctuations, and the embedding generates long- 
range correlations, that give to these large fluctuations 
their intermittent character. 



1.2.1 Scaling properties of discrete multiplicative cascades 

We symbolize the eddy cascade by cells, each cell being 
associated to a random variable Wi , called "weight" . All 
these random variables are assumed positive, independent, 
and identically distributed, with the following conserva- 
tive property: < W >= 1, where < . > denotes expecta- 
tion. The larger scale corresponds to a unique cell of size 
L = loXi, where 4> is a fixed scale and Ai > 1 is an adi- 
mcnsional scale ratio (very often for discrete models one 
takes Ai = 2). The model being discrete, the next scale in- 
volved corresponds to Ai cells, each of size L/Xi — foA"" 1 . 
This is iterated and at step p (1 < p < n) there are A^ cells, 
each of size L/\\ = i?rjA"~ p . There are n cascade steps, and 
at step n there are A™ cells, each of size £o, which is the 
smallest scale of the cascade (see Figure 1). To reach this 
scale, in agreement with Richardson's original idea d and 
also with more recent results on the localness of interac- 
tions in turbulence [8], all intermediate scales have been 
involved. Finally, at each point the energy flux writes as 
the product of n weights 9 : 

n 

e(x) = Y[ W p>x (2) 
P =i 

where W p . x is the random variable corresponding to po- 
sition x and level p in the cascade. Since each W PtX for 
different cells are assumed independent, their moment of 
order q > can be estimated as: 

n 

< «x)) q >= n < (^v) 9 >=< wq >n ( 3 ) 
i=i 

Introducing the total scale ratio 

A=^ = A? (4) 

between the outer and the inner scales, Eq. [5] gives fi- 
nally the moments of the small-scale flux, denoted here- 
after e\(x) to indicate its scale-ratio dependence 

< e\ >= X K(q) (5) 

where we introduced the cumulant generating function 
K(q) = log Ai < W q >. The conservative property < 
W >= 1 gives K(l) = and also < e >= 1. 

Equation [S] is a generic property of multifractal fields, 
obtained through multiplicative cascades. Depending on 
the model chosen, different forms of the pdf of W can 
be taken, leading to different analytical expressions for 
K(q). A log-normal pdf for W corresponds to a quadratic 
expression K(q) — ^(q 2 — q), with fi = K(2) = cr 2 /log2, 
where a 2 is the variance of log IF (this is detailed below). 
It can be also noticed that K(q) is (up to a logAi factor) 
the second Laplace characteristic function of the random 
variable log W (see the Appendix A) , showing that it is a 
convex function (see [Tfj]1. 

Let us mention here a specific property of multiplica- 
tive cascades: up to now, we have considered a cascad- 
ing process from large to small scales, called "bare" field 
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Fig. 2. The paths leading to 2 points separated by a distance 
r: they have a common path before a bifurcation. 



by Mandelbrot corresponding to the energy flux in 
turbulence. When the smallest scale is reached, the Kol- 
mogorov dissipation scale 77, the flux becomes the dissipa- 
tion. An interesting question is here to characterize the 
"dressed" field, corresponding to integrated small-scale 
dissipation as proposed by Kolmogorov and Obukhov 0] . 
Dressed quantities (e.g. observables) are a measure, up 
to a resolution scale r ^> 77, of the small-scale cascade, 
developed down to scale rj. In fact, the dressing may in- 
troduce divergence of moments |lllll8j but it is known 
that for moments below this critical moment qp, the scal- 
ing exponents of bare and dressed fields are identical. An 
interesting mathematical question is also the convergence 
of the measure when the total scale ratio A = L/rj — > 00 
|12|. But this point is not the main topic of our paper. 
Indeed, the dressing is usually seen as performed by the 
observer; we therefore consider here our stochastic process 
(presented below) as representative of the cascade devel- 
oped from a large to a smaller scale, i.e. a bare cascade. 
The question of the interrelation of dressed-bare fields and 
their respective statistical moments is therefore the same 
for our stochastic process, as for all multifractal cascade 
models (discrete or continuous). 



1.2.2 Correlation properties 

We now consider the long-range correlations generated by 
this cascade model. As originally shown by Yaglom [5], 
the correlation < e x (x)e x (x + r) > can be decomposed as 
a product of random variables, introducing their common 
ancestor as follows. The distance r corresponds in average 
to m steps, with A™ « r, so that the n — m first steps are 
common to the two random variables, whereas there is a 
"bifurcation" at the step m, and the two paths separate 
(see Figure 2). In the product of the two random variables, 
there are thus n — m identical variables and m different 
and independent. This gives: 

n n 

< e x (x)e x {x + r) >=< J[ W p , x J| W p >, x+r > 

p—1 p'—l 



n < w l > n < w v<* > 

p—l p=n —m+1 



X [] < Wp' lX+r > 
p f =n— 

< w 2 >"-™< w > 2m 



(6) 



Finally, introducing K(q) and recalling that A™ = A and 



A', 



r this gives: 



<e x (x)e x (x + r) >«A Jf ( 2 )r^ 



(7) 



Since A is fixed, this relation corresponds to a long-range 
power-law correlation with exponent /1 = K(2). Let us 
note that the condition A™ ~ r imposes that the scaling 
relation Q is exact for discrete scales only. This is an 
important limitation of discrete cascades, together with 
the fact that they are not strictly stationary. Relation {7J| 
can be generalized to provide, following the same line, the 
two-points correlations of moments of order p > and 

g >o nana; 

< e x (x)Pe x (x + r)« >« \K(P+<i) r -{K(p+<i)-K(3P)-K{<i)) 

(8) 

Since K(q) is non- linear and convex for multifractal dis- 
tributions, the exponent K(p + q) — K(p) — K{q) > 
quantifies the long-range power law correlations of multi- 
fractal measures. For p = q — 1, this yields the fi — K(2) 
exponent given above. 

On the other hand, we may also consider the corre- 
lation properties of the logarithm of the process (called 
generator), which possesses also interesting characteris- 
tics. Let us define the generators 7 and g of respectively 
the cascade process e and the weight W: 



7 = log e 
g = log W 

We consider the autocorrelation function of 7: 

C 7 (r) =< 7(21)7(0: + r) > 



(9) 



(10) 



As before, the path leading to 7(1) and 7(2; + r) may be 
split: for the last m steps the paths are different, whereas 
before this bifurcation, the paths are common. The main 
difference with the calculation leading to Eq.(jZJ) is that 
in Eq. (|10fl 7 is not the product but the sum of random 
variables. A simple calculation shows that C 7 (r) involves 
n 2 — n + m terms of the form < g > 2 and n — m terms of 



the form < g >. Let us note G 
a 2 =< (g - G) 2 >=< g 2 > -G 2 . 



=< g >=< \ogW > and 
Wc then have: 



C 7 (r) 



< (/ > 2 (n 2 — n + m)+ < g 2 > (n — m) 



ri2 2 
G n 



C- 



+ a (n— Ti 
logr 



logAi 



(11) 



where we introduced the constant C = G 2 n 2 + a 2 n. This 
corresponds to a logarithmic decay of the autocorrelation 
function of the generator |lfilll7| . Its Fourier transform 
gives the power spectrum of the singularity process j(x): 



E-y(k) 



(12) 
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corresponding to an exactly 1// noise |18U19U2()j . Prop- 
erties (0, (jHJ, ifTTJ) and (fl2")l may be used to check the 
multifractality of a stochastic process. Analogous proper- 
ties are verified for dressed multifractals and can hence be 
used to test numerical simulations or experimental data 
(for a link between log-correlations and dressed scaling 
properties, see also |27)ll2"T] V 

1.2.3 The log-normal case 

Up to now, the analytical expression taken by K(q) is only 
loosely constrained a priori: a convex function verifying 
if (0) = and K(l) — 0. We explicit here the generic 
example of multifractal cascade, given by the log-normal 
case, where generators are Gaussian. We choose here for 
simplicity Ai = 2. We take W — e 9 , where g is a Gaussian 
random variable of mean G and variance a 2 : g = en?o + 
G, where go is a centered and unitary Gaussian random 
variable. We have: 

< W q >= exp QcrV + (13) 

The conditions < W >= 1 and = K{2) = log 2 < W 2 > 
provide the constants, and finally the weights of the form: 

W = exp (V/ik^o - | log 2) (14) 

generate a discrete log-normal multifractal field with mo- 
ment function 

K(q) = ^(q 2 - q ) (15) 

This corresponds to log-normal statistics for e, with a 
variance of log e obeying Kolmogorov's initial hypothesis, 
given by Eq.QJ. Formula (|14f> can be directly used for 
discrete log-normal simulations. 

2 A causal log-normal stochastic evolution 
equation and its time and scale properties 

2.1 Introduction 

Up to now, a discrete cascade procedure has been pre- 
sented, and its characteristic properties given. In this sec- 
tion we are interested in the scale densification, introduc- 
ing more and more cascade steps, while maintaining con- 
stant the total scale ratio. We first recall below how this 
leads to log-infinitely divisible (log-ID) probability distri- 
butions. This considerably restricts possible pdf for con- 
tinuous scale multifractal cascades. On the other hand, it 
is also very useful to have a stochastic equation generating 
continuous scale multifractal field or process at our dis- 
posal. Such a stochastic equation would depend on some 
parameters and a scale variable for fields, or on a time 
variable for processes. A few recent papers have been de- 
voted to such topic, and we briefly present them below, 
before providing a stochastic equation corresponding to 
a causal log-normal multifractal process. We then study 



this stochastic equation: we show that this produces in- 
deed a multifractal log-normal process with the desired 
correlation properties, and we discuss further its differen- 
tial form, when differentiated versus time and versus scale 
variables. 



2.1.1 Continuous cascades lead to log-infinitely divisible 
processes and stochastic measures 

To densify the cascade described above, we keep the to- 
tal scale ratio A large but fixed; the continuous limit can 
be obtained by increasing the total step number n, hence 
Ai = A 1 /™ — > 1 + (see 19,22,23 ,24 ). It is then well-known 
that, in this limit, loge is an ID random variable (see 
Appendix A and 1 for ID random variables): continu- 
ous cascade models are log- ID [22 ,25,26, 27 . This restricts 
the eligible cascade models, since ID laws define a specific 
family of probability distributions. 

In order to provide stochastic equations and processes 
generating log-ID continuous multifractals, we have re- 
cently introduced ID stochastic integrals based on ID mea- 
sures |2"%ll2l?] (see also |3T)ll3*T] V ID measures are set func- 
tions, also called in the mathematical literature indepen- 
dently scattered random measures • This provided gen- 
eral expressions for the whole log-ID family. A more re- 
stricted but still quite general family is composed of log- 
stable processes, for which generating equations take sim- 
pler forms |18II14U29| . Among these, we consider here a 
special but important case: log-normal processes. We study 
this here because the stochastic integral involved is the 
classical Gaussian stochastic integral, whose statistical and 
differential properties are well-known; furthermore, log- 
normal multifractals are classical (and historical) exam- 
ples of multifractals. The results given here for this type 
of continuous cascade may be generalized later for other 
ID processes. 

Let us underline that the bare/dressed question men- 
tioned above for discrete cascades is still valid for continu- 
ous cascades; the convergence properties of dressed log-ID 
measures has been studied in [3TJIl3"3"| . 

2.1.2 The stochastic equation 

Until now, cascades have been presented in space; for the 
remaining of the paper, we consider cascades developing 
in time. The cascade is then seen as a stochastic process 
possessing long-range correlations for time increments r. 
The stochastic equation we study here is the following: 

e x (t) = A^ /2 exp (a 1 ' 2 £ (t + 1 - u)- 1/2 dB(u)j 

(16) 

where c > is a constant, A ^ 1 is the total scale ratio, 
t is time and B(t) is a Brownian motion. This was al- 
ready given in |29| and in other less explicit forms in |18l 
rTMrl] but its statistical and differential properties have 
not been studied in details in these papers. This equation 
corresponds to the exponential of a Gaussian stochastic 
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integral, this integral being analogous to a fractional inte- 
gration of a Gaussian process, i.e. a fractional Brownian 
motion |15j . The only difference with a formal fractional 
integration comes from the fact that the domain of inte- 
gration is finite (for large values of A it is a wide domain, 
but still finite). The exponent of this fractional Brownian 
motion is not a free parameter. Let us recall that for a 
self-similar parameter H the kernel in the integration for 
a fractional Brownian motion is on the form (t — u) 11 ^ 1 ^ 2 , 
so that the kernel introduced above corresponds to a self- 
similar parameter H — 0. This is thus a very special case 
of fractional Brownian motion. This equation is also causal 
(the past does not depend on the future) and can be used 
to generate sequentially a log-normal process with no up- 
per limit in the number of points. We study below in de- 
tails the properties of this equation. 

It is first easily verified that this process is stationary: 



(17) 



where = means equality in distribution. The statistical 
moments are then independent of time, and are estimated 
below: 



/m\o \-au/2 W 1/2 P (t+l-u)- 1 ' 2 dB(u) 



> 



exp 



\xq 2 f du 



t+ 



!_ A t + 1 - U 



A-^A^ 



(18) 



giving the requested scale-invariant properties, Eq.(|SJ) and 
Eq. Ijl5|l (see Appendix B for the second characteristic 
function of a stable stochastic integral). We explore be- 
low other properties of this equation in time and in scale- 
space. 



2.2 Time statistics of the log-normal stochastic 
process 

2.2.1 Correlations 

We consider here the generalized correlation function 

C tx (p,q,r) =<exWe X (t + T)* > (19) 

for moments of order p > 0, q > and for 1 < r < A. 
The stochastic integrals are split in order to separate the 
overlapping integration domains, corresponding to inde- 
pendent random integrands. This gives: 

C £ >,g,r) = A^^ x 

S+l+T-A 

x < cxppy/JI J (t+1- u)- 1/2 dB(u) > 

t+l-X 

t + T 

x < exp qy/JI J (t + l + r- uy l/2 dB(u) > 



Z 

x<exp^7i J (p(t + 1 - u y 1/2 + 

t+l+r-A 

+q(t + 1 + T- u)- 1/2 )dB(u) > 



t + l + T-X 

2 '' du 



t+T 



x exp 



t+l-A 



du 



g \i 

2 J t+l+T-U 

t 



= i lA ^+«) (_A_) V ( T + i)4- 

where 1\ is the last integral to evaluate: 



(20) 



Ji = exp 



{pit + l-u)- 1 ' 2 + 



t+l+T-X 

+q(t + 1 + r - u)- 1 / 2 ) 2 




x exp 



du 



gV 

2 J I + "I + T II 

t + l + T-X 



h(X-r)- 



T+ 1 



(21) 



with finally 



h = exp fipq 



du 



t+l+T-A 



y/(t+l-1l)(t+l+T-u) 



2fipq 



\ 1 + V7TT J 

where we have used the identity 



(22) 



dx 



\/x{x + a) 
Eqs.JUJ, (T5TJ) and HU give: 

C ex (p, q ,r)=\ K W+ K M[ 



= 2 log (■sfx + Vx + a) (23) 



k(p,q) 



i + VtTT 



(24) 



with k(p, q) = 2 (K(p + q) - K(p) - K(q)). Whenever 1 < 
t <C A this gives finally: 

C ex (p,q 7 T) « xK(P+i) T K( P )+K( q )-K( P + q) (25) 

having the same asymptotic form as the generalized cor- 
relations for discrete cascades (Eq. (JHJ). 
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On the other hand, we may introduce the logarithm, 
called here singularity process 

7x(t) =bge A (t) = -f logA + 
t 

+ A t 1 < /2 J (t+l-u)- 1/2 dB(u) (26) 

t+l-A 

and estimate its autocorrelation function for 1 < r < A — 1 
(using Appendix B): 

C 7A (r) =<Jx(t)7x(t + r) >= (§log'A) 2 + 

+M< J (t + l-u)- 1/2 dB(u) x 
t+l-A 

t+T 

{t+l + T-u)- 1/2 dB(u)> (27) 

t+l-A+r 

As before, the stochastic integral can be split into two non- 
overlapping domains. Then, using Eq. l|59|) from Appendix 
B, we have: 



Gt x (t)= £ log A 



■ fi j (t + i- u y 1/2 x 

t+l—X+T 

x(t + l + r - U y l/2 du (28) 
and using again Ea.l|2*3"]l. this gives: 



C^{t)= (flog A) + 2 M log^ 

~ A\ — [l log T 



(29) 



where A\ is a constant depending weakly on A. The last 
line used as before the assumption 1 <C r <C A. Equa- 
tion (|29J) is very close to Eq. Qlljl which was obtained for 
discrete cascades. 

The stochastic process introduced in Eq. I|16[l has thus 
all the characteristic properties of multifractal cascades: 
multiscaling of the 1-point statistics, together with long- 
range power-law correlations for the process and log-corre- 
lations for the logarithm of the process. We consider below 
other properties of this equation: composition rules and 
differential forms. 



2.2.2 Differential properties 

Let us write stochastic equations for the processes e\(t) 
and 7a (t). First, let us consider the process Vi{t) = f {t + 
l—u)~ l ^ 2 dB(u). Applying the result given in Appendix B, 
V(t) is a martingale process (with respect to the natural 
filtration) having the differential form: 



The same formula applied to the process Vi (t) = J Q (t + 
A — u) _1 / 2 d.B(u) gives 



dV 2 (t) = \- 1 / 2 dB{t) + W 2 (t)dt 

w 2 (t) = ~y*(t + \-u)- 3 / 2 dB(u) 



These relations give then the stochastic process 7a (i), hav- 
ing as differential form dy\(t) = /i 1/2 (dV"i(t) - dV 2 (t+ 1 - 
A)): 

d lx {t) = li 1 ' 2 (dB(t) - \-^ 2 dB(t + 1 - A)) - \W{t)dt 
W{t)^^/ 2 J (t + l-u)-^ 2 dB(u) 

t+l-A 

(32) 

For large scale ratios (A> 1), the term A 1 l 2 dB(t + 1 — A) 
may be neglected, since it is much smaller than dB(t), so 
that the differential takes a simple form. 

We consider now the stochastic differential equation 
for e\(i) = exp7A(i). Ito stochastic calculus gives (recall 
that dX(t) = X{t)dY{t) + \X(t)dt for X(t) = expF(i)): 

de x (t) = fi 1/2 e x (t) (dB{t) + i(l - W(t))dtj (33) 

where we have still neglected the term \~ x / 2 dB{t + 1 — A), 
and W(t) is given in Ea. f32p . This provides the stochastic 
differential equation having as a solution the multifractal 
lognormal process studied here. It corresponds to a multi- 
plicative process with a stochastic drift term. This may be 
compared to other families of stochatic differential equa- 
tions; it can also be used for theoretical analyses, or to 
estimate (theoretically or numerically) the predictability 
properties of such processes. Indeed the long-range cor- 
relation properties of such process correspond to a long 
memory that can be exploited to optimize its predictabil- 
ity. This was already done in another context |14| and 
could be pushed further for the present stochastic pro- 
cess. We do not develop further this point here, and let 
such analysis for future studies. 



2.3 Scale statistics of the log-normal stochastic 
process 

In the previous section, we have considered the time prop- 
erties of the stochastic evolution equation, for a given (and 
large) value of the scale ratio A: time correlations and time 
stochastic differential equation. In this section we consider 
scale properties, assuming a fixed time t and considering 
the development of a cascade at a fixed point, when A 
increases. We then define 



E t {\) = e x (t) 



(34) 



( dVi(t) = dB{t) + Wi(t)dt 

\W l (t) = - 1 1 f (t + l-u)-^ 2 dB{u) 



(30) 



We take A as a continuous variable; in some instances we 
consider the log of the scale ratio and take as scale variable 
P = log A. 
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2.3.1 Semigroup composition properties 

Let us first recall that, in consequence of their construc- 
tion, discrete cascade, when considered at a fixed position, 
can be considered as a discrete Markov process in scale 
(see e. g. EU): since successive multiplicative weights are 
independent, the product of n + m weights can be decom- 
posed into the product of n weights times the product of 
m weights, possessing a semi-group property. 

We consider here a semigroup property hidden in our 
stochastic evolution equation. 

E t {\iX 2 ) = (AiA 2 )-^ 2 x 

t+l-Ai 

- cxp I p 1/2 [ (t + 1- u)~ 1/2 dB(u) 



t+l-AiA 2 



exp L 1 / 2 J (t + l-u)- 1/2 dB(u) J 

\ t+l-Ai / 



x exp 



/ t/Ai 



,V2 



\ t/Ai-A 2 + 



Ai 1/2 {t + X^l ~ v)y 1/2 dB(v) 



^ t (Ai)ex P U 1/2 / (*' + l-v)-V 2 dB(v) 

V t' + l-A 2 

Et{\i)E* (A 2 ) 



(35) 



In the first equation the integration domain has been split 
in two non-overlapping domains, so that the corresponding- 
random variables are independent. In the second equation, 
we have introduced a change of variable (v — 1 = 3 y^ : ) 

and the fact that dB(Xv) = X 1 ' 2 dB(v). The third and 
fourth equations are simple factorizations (with t' = 
using the definition of E t {X). As a result, Eq.J2SJ) is a 
decomposition property "in distribution" , for a fixed time. 

Let us note that usually in papers studying Markov 
properties of the cascade in scale, the time or space vari- 
able do not appear. This is justified by the fact that the 
process is stationary, and that one often considers the 
probability instead of the process itself. To see this, let 
us introduce tt p {x) — Yx(\ogE t (e p ) — x). Then taking the 
log of Ea. (|35[l . one obtains an equality in distribution for 
the sum of random variables, leading to a convolution of 
their probabilities: 



'P1+P2 



'Pi 



'P2 



(36) 



where pi = log A, . 

This recovers the convolution semi-group property whi- 
ch can be found also in several recent papers |36U27U35| . 
In these papers, the convolution property has been shown 
to be directly linked to a decomposition property (origi- 
nally introduced in [37]) of the pdf of the cascade process. 
This can also be obtained considering that the process 



f(p) — \ogE t (e p ) is a continuous Markov process with ho- 
mogeneous transition probability of the form Q r (x,y) = 
P r (/(/° + r) — y\f(p) = x) = q r {y/x). Then Chapman- 
Kolmogorov equation for transition probabilities can be 
written on the following form 



q r+s {x) = I q r (-J q s (a)da; a > 



(37) 



This equation is a general property of log-infinitely divis- 
ible (log-ID) processes, but since characteristic functions 
have simple expressions for such processes, it is often more 
convenient (as we do in the present paper) to characterize 
them using characteristic functions instead of probability 
distributions. 

In a series of recent papers, Friedrich, Peinke and col- 
laborators 40 41 42 have experimentally shown that, for 
a fixed position, the velocity or dissipation statistics in 
turbulence are Markov processes. Their Fokker-Planck ap- 
proach can be shown to be linked (in case of lognormal cas- 
cades) to Castaing's semi-group decomposition approach 
38,39,36 . Below we consider in this framework different 
properties in scale, leading to a Langevin equation. 



2.3.2 Differential properties and a Langevin equation 

The series of papers about experimental evidence of a 
Markov process for the energy cascade in turbulence have 
considered either the velocity field or the dissipation field, 
and focused on a Fokker-Planck or Langevin (differen- 
tial) equations (for a link between log-normal multifractal 
statistics and Fokker-Planck approach, see |HE])- Here we 
give the theoretical expression obtained for our process 
and compare our result with experimental results pub- 
lished for the dissipation field in turbulence. 

A stochastic (Langevin) equation has been obtained 
by Marcq and Naert 03| f° r the logarithm of the cascade 
process (more precisely the centered variable), involving 
the logarithmic scale ratio variable p = log A. We check 
here what this gives for our process. First for its generator, 
we have: 



lt (p)=\oge x (t) = -1p + p 1 / 2 G(eP) 
G(x)= J (t+l-uj-^dflfu) 

t+l-x 

This is easily differentiated and provides: 

1dp + p 1 / 2 ePdG{eP) 



\ dG(x) 



-1/2 



dB(t +l-x) 



(38) 



(39) 



This may be more simply written introducing another 
Brownian motion Bo(X) = —B{t + 1 — A) (recall that the 
time is frozen here) , and writing the equation for the cen- 
tered singularity 70 = 7*— < 74 >: 



djo(p) = VJI^dBo{e p ) 



(40) 



where we used the property < G(x) >— and hence < 
7t(p) >= — 77 p. For the process Et(X) itself, this gives the 
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following equation: 

dE t {\) = E t {\) (^-{1 + ^\V* d B (\fj (41) 

This gives the Langevin equation in scale for our contin- 
uous multifractal process (for a fixed time). 

Our results can be compared to previous studies, some 
of them done in purely experimental grounds. Let us first 
note that in some cases (see 0J) the study is done con- 
sidering a Fokkcr-Planck equation with drift and diffusion 
coefficients D\ and D 2 . The associated stochastic differ- 
ential equation (Langevin form) is: 



dj(p) = D 1 (p)d P + ^/2D 2 (p)dB(p) 



(42) 



so that our results given in Ea. H4U|) for the centered vari- 
able provide an equation analogous to Eq. 142f) with: 



D 1 (p) = 
D 2 (p) = \h & p 



(43) 



Our stochastic equation is thus analogous to a zero drift 
and and scale-dependent diffusion coefficient, but with a 
subordinated Brownian motion B${e p ), instead of the pro- 
cess itself. Let us underline that this applies to a bare 
cascade, developing from large to small scales, at a fixed 
time. This can be compared to the results given by Cleve 
and Greiner for a discrete bare cascade, where they 
also find a zero drift coefficient. On the other hand, ex- 
perimental results on the Markov property of the cascade 
correspond to properties of the dressed cascade and there- 
fore may not be related directly to properties of the bare 
cascade. This difference was already underlined in 
where it was shown numerically (from a discrete cascade 
simulation) that the dressing introduces a non-zero drift 
term. The theoretical reasons to obtain a Markov prop- 
erty for dressed quantities is also less straightforward than 
for bare ones. Nevertheless, we may compare the values 
given by our bare model with dressed experimental results. 
Naert et al. |41) have experimentally obtained a non-zero 
drift value and a constant D 2l later Marcq and Naert 03| 
have found a scale-dependent diffusion term of the form 
D 2 (p) = D e 25p with D « 0.01 and S « 0.4. This is 
closer to our result: Do is clearly too small in 02| com- 
pared to our expression (recall that experimentally, p is 
close to 0.25, leading to Dq w 0.12) but the p-dependence 
is almost the same since we would have S = 1/2. In any 
case we do not expect here to recover the same values, 
since (i) there may be a bare-dressed difference in the co- 
efficients and (ii) the discrepancy between our theoretical 
expression and experimental results are also likely to come 
from the fact that we took for simplicity (as a first step) a 
log-normal cascade here, whereas the real picture in tur- 
bulence seems to be more involved (see e.g. |H]L 

Let us underline once more that this scale description 
is valid only for a fixed position in time or in space, since 
the long-range correlations of turbulence in time or in 
space are not compatible with a Markov process. A more 
complete representation of the cascade process is then - 



as we do in the present paper - to provide a time-scale 
description. In this time-scale space, the process involved 
is more complex than being simply Markovian. 



2.3.3 Interpretation of multiplicative cascades as stochastic 
continuous products 

Discrete cascades correspond to a finite product of ran- 
dom variables. For convenience reasons, the densification 
has been interpreted above as the exponential of a stochas- 
tic integral. In fact, the natural way to densify the cascade 
would be to introduce more and more intermediate scales 
in the finite product, while keeping the total scale ratio 
fixed. This corresponds to a continuous product (for de- 
terministic continuous products see |45j ) . Stochastic con- 
tinuous products (also sometimes called stochastic multi- 
plicative integrals) have been introduced by McKean , 
and developed by e.g. Ibero 02], Doleans-Dade 0B] and 
Emery |49j in the framework of continuous martingales. 
As first argued in [2E], continuous products are necessary 
when one needs to introduce zeroes (such as in a contin- 
uous generalization of the discrete /3-model HOD; or when 
generalizing the ID approach to matrix or tensorial cas- 
cades. Stochastic continuous products provide in fact an 
interesting mathematical framework to study continuous 
cascades, since some results obtained in this field may be 
applied to multifractal cascade processes. Some first re- 
sults in this approach are given below. 

Among the results obtained by Ibero 02] is the fol- 
lowing: under general conditions on the square matrices 
functions / and e (they can be stochastic functions), the 
following multiplicative stochastic integral (denoted as p|, 
see Appendix C for a definition) is converging 

t 

z(t) = P| exp (f(s)ds + e(s)dB{s)) (44) 
o 

and is the solution of the following stochastic differential 
equation: 



1 



dz{t) = z{t) (/(t) + -e 2 (t))dt + e{t)dB{t) 



(45) 



in the algebra of square matrices of order n (B is a matrix 
Brownian motion). This framework is useful for cascades 
of matrices, which can be used to consider vectorial multi- 
fractal processes. In this line, the generalization of Ea. (|16(l . 
is, introducing continuous product: 

t 

e x (t) = \- p/2 f] e^(p 1/2 {t + \-u)- 1 ' 2 dB{u) S j 

t+l-A 

(46) 

which will give an equation analogous to Eq. 1)33(1 for 
square matrices. Indeed, for ID processes, this equation 
gives back Eq. I|16|) , while for the general multi-dimensional 
case where B is a matrix Brownian motion, this stochastic 
continuous product can be shown to correspond to a ma- 
trix multifractal process. This means that each element of 
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such matrix is a ID multifractal process, while they also 
possess two-by-two scaling properties. We do not explore 
further here the multifractal and correlation properties of 
Eq. H46J) : this will be done elsewhere. Neverthelsess, it is 
clear that it would be interesting to develop such theo- 
retical framework in order to perform data analysis, such 
as for example for analyzing multidimensional data series. 
We discuss here a potential application of such approach: 
these cascades of matrices may be used to model intermit- 
tent anisotropic tensors in turbulence, such as the dissi- 
pation rate £ij(t) where it is defined as 



2v < 



dui duj 
dxk dxk 



> 



(47) 



Indeed, according to Kolmogorov's 1941 hypothesis, for 
high Reynolds number flows the cascade process is as- 
sumed to "wash" the details of the large-scale flow so that 
after enough cascade steps the dissipation tensor becomes 
isotropic. In fact experimental results in homogeneous tur- 
bulence indicate that this may not be true |51) . Several 
deterministic closure models have been developed to take 
into account this anisotropy for homogeneous and non- 
homogeneous turbulence (see e.g. |52U53| ). This type of 
model are fully tensorial and can be numerically solved to 
provide the prediction of some statistical quantities (such 
as mean and second moments of the velocity field). But 
they need too drastic hypothesis (in particular they are 
not compatible with the long-range correlation properties 
of turbulence) and they provide, at best, adequate predic- 
tions for moments up to order 2 (see a general presenta- 
tion on this in jSU)- On the other hand, turbulent cascade 
models are able to reproduce experimental fluctuations of 
the velocity field up to moments of order 7 (see but 
up to now they can deal only with ID processes. Because 
of this limitation, they have offered until now very limited 
industrial applications. If turbulent cascade models could 
be adapted to vectorial or tensorial frameworks, it could, 
in a relatively short time, lead to much more reliable pre- 
dictions for industrial flows. This promising research di- 
rection is not developed further here. 



3 Simulations 

We perform here numerical simulations of discrete and 
continuous cascades. The aim of these simulations is to 
check the theoretical relation given in previous sections, 
and to compare discrete and continuous cases. 



3.1 Simulation of discrete log-normal cascades 

We perform here a numerical simulation of a discrete log- 
normal cascade with Ai = 2 and weights given by Eq. l|14l) . 
This is illustrated here with n = 12 cascade steps, corre- 
sponding to a total scale ratio of A = 2 12 = 4096, and a 
value of the intcrmittency parameter close to experimen- 
tal results for the dissipation: \x = 0.25. Figure 3 shows 



E 24 




500 1000 1500 2000 2500 3000 3500 4000 



X 



Fig. 3. A sample of a log-normal discrete cascade with total 
scale ratio 2 12 . The generator 7 (a correlated Gaussian noise) 
and the multifractal field e are shown. 




Fig. 4. The autocorrelation function of the generator, esti- 
mated from 100 realizations of a discrete log-normal cascade 
with n = 0.25. As expected this shows a logarithmic decrease 
over the whole scaling range. Dotted line: predicted logarithmic 
decay of the form — fj, log T. 



the resulting field for the cascade process e and the sin- 
gularity 7. It is visible that e is highly intermittent while 
7 is a correlated noise. This correlation is illustrated in 
Fig. 4, showing C 7 (r) versus logr: the logarithmic decay 
is extremely well verified for almost 3 orders of magnitude, 
with a slope of about 0.24, close to the predicted value, 
which is a 1 1 log 2 = fi. Figure 5 represents the power spec- 
trum of 7 and of e in log-log plot. These are both scaling, 
with respectively theoretical slopes of -1 and —1 + [i: see 
Eq. (|12J) and the Fourier transform of Eq. 10. Experimen- 
tal results are very close to theoretical power laws for the 
whole range of scales considered. 

Figure 6 is a direct test of the scaling of the result- 
ing field, for various order of moments: 100 realizations 
(each of length 4096) have been taken into account for the 
estimation of statistical moments. The dressed moments 
are represented in log-log plot, for different values of q; 
the straight lines confirm a scaling law of the form given 
by Eq03for the dressed field. The slope of these straight 
lines give experimental estimates of K{q) 1 and the result- 
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Fig. 5. The Fourier power spectrum of 100 realizations of a 
discrete log-normal multifractal field and of its generator. As 
expected, over the whole scaling range, the generator displays 
an exactly / power spectrum, and the multifractal field a 
Dotted lines provide theoretical behaviours. 




1 10 100 1000 10 4 A, 

Fig. 6. The scaling of the moments of order 0.5, 1.5, 2, 2.5, 3, 
3.5, and 4 (from below to above), estimated for 100 realizations 
of the discrete log-normal multifractal cascade. The different 
straight lines indicate the accuracy of the scaling property. The 
slopes of these straight lines give estimates of K(q). 



ing curve is shown in Fig. 7, compared to the theoretical 
bare expression § (q 2 - q). There is an excellent agree- 
ment until a critical moment q s ~ 2.4, above which the 
estimated function is linear: this corresponds to a maxi- 
mum order of moment that can be estimated with a finite 
number of samples, since there are not enough realizations 
for an accurate estimation of the scaling exponent corre- 
sponding to larger order of moments. This may also corre- 
spond to the critical moment above which moments of the 
dressed field diverge; in this case experimental estimates 
are not infinite, but the exponent K(q) becomes linear 
(see |5(:i| for a discussion on bare/dressed exponents in an 
experimental context). In any case, experimental analy- 
sis shows that the simulated field has scaling moments, 
and as expected for low orders of moments, the moment 
function is close to the theoretical bare curve. 



Fig. 7. The K(q) function estimated from 100 realizations of 
the discrete multifractal simulation (open squares), from 106 
points of a continuous simulation (filled dots), compared to the 
theoretical curve given by Equ. (17) (dotted line). There is an 
excellent agreement until a critical moment corresponding to 
the maximum moment that can be accurately estimated with 
a finite sampling. 




1000 2000 3000 4000 5000 t 



Fig. 8. (a) A sample of the Gaussian noise and of the gen- 
erator; the generator is a moving average over the past of the 
noise, and corresponds to a correlated Gaussian noise, (b) A 
sample of the multifractal process, corresponding to the expo- 
nential of the generator of Fig. (8a) . 

3.2 Simulation of continuous and causal log-normal 
processes 

Here we provide a stochastic simulation of the process 
given by Eq. (jTB)l . choosing A = 5,000 and fi = 0.25. 
This is done sequentially, generating the noise and a mov- 
ing average over this noise. An arbitrary large number of 
points can be generated with this equation. We represent 
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C(T) 




-0.5 t 1 

1 10 100 1000 

T 

Fig. 9. The autocorrelation function of the generator, esti- 
mated from a continuous stochastic process with ii = 0.25. As 
expected this shows a logarithmic decrease for large time lags. 
Dotted line: predicted logarithmic decay of the form —fx log T 
with an experimental [i value of 0.23. 



E(f) 




1 10 100 1000 f 



Fig. 10. The Fourier power spectrum of the continuous 
stochastic process. As expected, over the whole scaling range 
the multifractal field displays a regime with an experi- 

mental value of fi — 0.23 (dotted line). 



in Fig. 8a-b a sample path of 5,000 points for the process 
and its generator: Fig. 8a shows the Gaussian uncorrelated 
noise, together with the generator process, which is a mov- 
ing average over this noise (with a specific kernel). This 
correlated noise has the same visual aspect as the gener- 
ator obtained in Fig. 3 for discrete cascades. Figure 8b 
shows the field e itself (the exponential of the generator): 
it presents as expected a huge intermittency. 

To estimate the statistics of the path, we have gener- 
ated 10 6 points. The scaling of their power spectra, to- 
gether with their correlation properties, have been veri- 
fied. The logarithmic decay of the autocorrelation function 
of the singularity process is shown in Fig. 9: it is nicely 
verified with a parameter of \i = 0.23 close to the theoret- 
ical value. For short time lags (smaller than about 10) the 
experimental curve leaves the straight dotted line corre- 



sponding to the fit; this is expected since the logarithmic 
decay was shown to be an asymptotic law for large enough 
time lags. Figure 10 shows the power spectrum of the field 
e which is scaling with a fit of the form f^ 1 as expected, 
with an exponent of —0.77, close to the expected value. 
The resulting moment function K(q) is also shown in Fig. 
7, as a comparison with the theoretical expression and 
the function estimated for the discrete simulation. The 
value corresponding to the continuous simulation is very 
close to previous results for the discrete simulation, and 
the same comments about bare/dressed scaling exponents 
apply. These different results confirm the different asymp- 
totic expressions given for 1 <C r <C A and the multifractal 
properties of the stochastic simulation. 



4 Concluding comments 

We have explored here the statistical properties of time se- 
ries generated by a continuous and causal stochastic equa- 
tion taking the form of the exponential of a stochastic in- 
tegral. We have first shown that such generating equation 
indeed corresponds to a log-normal multifractal: this was 
verified with one point and two points statistics, either an- 
alytically and also numerically, in the asymptotic regime 
(for large enough time lags) . We have also shown the com- 
position rule of such equation, and considered its differ- 
ential properties. This equation depends on 2 parameters: 
the intermittency parameter /j, and the total scale ratio 
A. For a fixed time and a variable scale ratio parameter, 
another type of process is obtained; the differential with 
respect to scale parameter leads to a Langevin equation 
for the logarithm of the process. We have linked this ap- 
proach to experimental results about Langevin equations 
for the dissipation in turbulence. Of course, this does not 
indicate that log-normal cascades are the best models for 
turbulent intermittency. This is rather the exploration of 
an important and classical multifractal model; other con- 
tinuous multifractals are associated to more complex gen- 
erating equations having less straightforward differential 
properties. 

This equation for log-normal multifractals is neverthe- 
less believed to be of great interest since a quite sim- 
ple stochastic equation generates a causal and continu- 
ous multifractal process, that can be used as a rather re- 
alistic approximation for various applications. It is also 
the first explicit continuous equation of its kind. This ap- 
proach is finally believed to open the door to many in- 
teresting generalizations, including cascades of matrices 
using a stochastic continuous product framework, which 
is only quickly evoked here. 
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Appendix A. Characteristic functions and in- 
finitely divisible distributions 

We recall here some useful properties of characteristic 
functions, infinitely divisible distributions, Gaussian and 
Levy-stable random variables. Useful references for these 
topics are Feller rjOj , Janicki and Weron , and Samorod- 
nistky and Taqqu [55]. 

Let X be a random variable with probability density 
p(x) . The usual characteristic function is the Fourier char- 
acteristic function 4>{k) defined as 



<j>(k) =< e 



ikx 



> 



(48) 



Since when studying multiplicative cascades real moments 
are introduced, we use here Laplace characteristic func- 
tions, assuming that the positive tail of the probability 
density decreases fast enough for the integral to converge: 



/(,)=< 



(49) 



One thus have f(q) — <fi(—iq) where 4>{q) are given in the 
usual tables of characteristic functions. 

Infinite divisibility is a property that has no simple ex- 
pression for probability densities. Instead, this property is 
very simply characterized using characteristic functions. 
A probability distribution is said to be infinitely divisible 
(ID) if its characteristic function / has the following prop- 
erty: for any n integer, there exist a characteristic function 
/„ such that: 

/" = / (50) 

In such a case one can introduce the second characteristic 
function & — log /: 



V(q) = log < e qx > 



(51) 



Let us consider a second characteristic function ^(q). If for 
any a > 0, a^{q) is still a second characteristic function, 
then the distribution associated to ^(q) is ID. 

Let us give two examples: a Poisson distribution with 
mean a has a second characteristic function of 



V(q) = a(e q - 1) 



(52) 



Levy-stable random variables can be defined as follows: let 
{Xi)i—x.. n be independent random variables of the same 
law, and n any integer. Then the law is stable if there ex- 
ists a n > and b n such that a n J27=i — b n has the same 
law as the Xi. Strictly stable random variables are the ones 
for which b n — 0, corresponding to centered < Xi > = 
variables. Using the second characteristic function, it is 
easily seen that Levy random variables, for which 



one recovers gaussian random variables with a 2 — r 2 /2. 
Levy-stable random variables are characterized by the ba- 
sic parameter a, limited to the range < a < 2. One 
may note that Levy-stable variables can take high values 
with hyperbolic probability tails, so that for the Laplace 
characteristic function to converge, one must consider only 
skewed Levy variables, for which the hyperbolic tail cor- 
responds to negative values only. In this case, the Laplace 
characteristic function is defined for moments q > 0. 



Appendix B. Stable stochastic integrals 

A useful reference for this topic is e.g. Janicki and Weron 
|57j . The rescaling properties of stable random variables 
can be used to define stable stochastic integrals. Let us 
consider again Y = X)"=i^j where (Xi)i=i., n are inde- 
pendent strictly stable random variables having the same 
law, with parameters a and r x . Then one has &y{q) = 
nfy x {q) — nF x q a showing that: 



r v = n 1 ^ 



r x 



(54) 



This property can be used to build stable random mea- 
sures and in a consistent way, stable stochastic integrals. 

For an interval of width dx, M{dx) is defined as a 
strictly stable random measure, i.e. a strictly stable ran- 
dom variable of index a and of scale parameter Z~m = 
(dx) 1 /". It has the following second characteristic func- 
tion: 

*M(q) - (r M qT =q a dx (55) 

For a positive valued function F, such that F a {x)dx 

exists, the stochastic integral I = J F(x)M(dx) is a ran- 
dom variable defined as follows: it is the strictly stable 
random variable of scale parameter 



r, 



\ 1/a 

F a (x)dx 



(56) 



which can be rewritten in the following way, using its sec- 
ond Laplace characteristic function: 

log < e q Sl F ^ M ^ >= [ [ b F«( x )dx) q a (57) 



This fully characterizes stable stochastic integrals. In par- 
ticular, for a — 2, the Gaussian stochastic integral / = 

F{x)dB{x) is still a Gaussian random variable with the 
variance 

b 

crj=[ F 2 {x)dx (58) 



¥(q) = r a q a 



(53) Gaussian stochastic integrals have also the property: 



are stable with a n = n 1 /". The parameter r > is called 
the scale or dispersion parameter; it plays the same role 
as the variance for gaussian random variables: for a — 2, 



< 



F{x)dB(x) / G{x)dB{x) >= 



E-2 



F(x)G(x)dx 
(59) 
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Indeed, if E± D E^ = the two integrals and independent 
random variables and whenever E\ U Ei 7^ only the 
intersection contributes to the correlation. 

Let us finally consider the stochastic process defined 

as: 

W K (t) = f K(t,s)dB(s) (60) 
Jo 

where the kernel K possesses some regularity conditions: 

K(t,s) = 0; s>t (61) 

pT / \ 1/2 

J I J d t K(u,s) 2 ds) + 

+ [ K 2 (s,s)ds < 00 (62) 
Jo 

then Wk (t) is a scmimartingale with respect to the usual 
filtration and has the following representation 

dW K {t) = K{t, t)dB[t) + W dtK {t)dt (63) 

See for more details Carmona and Coutin |59| . 



Appendix C. Stochastic continuous product 

Stochastic multiplicative integrals or continuous products 
have been introduced in the field of continuous martin- 
gales (see McKean [JS], and developed by e.g. Ibero |4"7) . 
Doleans-Dade [JH| 1 Emery and also Karandikar |60|L 
Below we propose a personal and short didactic introduc- 
tion to such objects. 

Infinitely divisible distributions are naturally associ- 
ated to random measures, which can also be called ran- 
dom additive set function. The discrete addition of ID 
random variables, when dcnsificd, leads to stochastic in- 
tegrals. The same way, one wants here to density a discrete 
multiplication of random variables, to obtain a continuous 
product. This can be done introducing random multiplica- 
tive set functions (RMSF). We associate with a measure 
m(A) for a set A ^ 0, a multiplicative set function n(A) 
defined as n (A) = e m ' A '. Then for A and B such that 
A n B = 0, one has n(A U B) = n{A)n{B). A random 
multiplicative set function N can then be defined for a 
set A: N(A) = e M ' ' and is characterized by the con- 
trol measure m(A). It is clear that for A and B with a 
void intersection, one has multiplicative random variables: 
N(A U B) = N(A)N{B). 

This is useful to introduce a random continuous prod- 
uct. Let us take non- intersecting intervals Ai and intro- 
duce the random variable P: 

n 

P = Y[N(Ai) (64) 
i=i 

This can be generalized to a continuous limit, taking in- 
tervals of the form Ai = [a + i^, a + (i + 1)^], whose 



reunion is the interval [a, b] and whose width decreases 
with n. This way one has the limit: 

b n ~ 1 f h — h — \ 

f| N(dx) = lim TT N I [a + i -, a + (i + 1) -] 

i=0 v 7 

(65) 

One obtains a stochastic continuous product, as the limit 
of a product of random multiplicative set functions, noted 
P). It has the following basic property: 

< (f] A N(d x )Y >= e™^* ® (66) 

which can be used for continuous cascades (!^o (<?) is a ref- 
erence second Laplace characteristic function). 
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